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We discuss the importance of using partially quenched theories with three degenerate quarks for extrapolating 
to QCD, and present some relevant results from chiral perturbation theory. 



Simulations of lattice QCD have, to date, 
largely used the quenched (Q), or partially 
quenched (PQ), approximations. These approxi- 
mations introduce unphysical artifacts, some of 
which can be analyzed using chiral perturba- 
tion theory (ChPT) |l|-^. In particular, it has 
been found that the quenched theory has un- 
physical singularities in the chiral limit, and that 
these persist in the PQ theory, albeit in weak- 
ened form ■ These results serve as a warning 
against relying on ChPT when straying too far 
from unquenched theories in which valence and 
dynamical masses coincide. 

In this talk we discuss a more constructive as- 
pect of partially quenched QCD. We focus on 
PQ theories with three light dynamical quarks, 
where "light" means that the theory can be stud- 
ied using ChPT. Such theories have the same 
quark complement as QCD, although with dif- 
fering masses. When simulating such theories 
it is possible to push to much lighter valence 
quarks than dynamical quarks, and thus to map 
out a "rectangular" region in the space of theo- 
ries as sketched in Fig. |l]. Our aim is to make 
use of the extra information contained in the PQ 
data to learn about unquenched theories, and 
in particular about real QCD. The key point is 
that, as long as one stays far enough away from 
the my = axis, ChPT rehably predicts the 
functional dependence on valence and dynamical 
quark masses, and so provides a unified descrip- 
tion of both the PQ and fully unquenched regions 
of the parameter space. 

An observation central to this work is that the 
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Figure 1. Schematic representation of the space 
of PQ theories. 

unknown coefficients appearing in the chiral La- 
grangian for PQ theories depend only on the num- 
ber of light dynamical quarks, and not on their 
masses 



One way of viewing the extra in- 
put that PQ theories provide is that they allow a 
cleaner separation of the non-leading, 0{p*) co- 
efficient^. This is because some GL coefficients 
appear multiplied only by valence masses, others 
only by dynamical masses (see Eq. (^). 

The picture that emerges is that, by simulat- 
ing in a region of moderately light valence and 
dynamical quarks, and varying independently the 
masses of both, one can determine all the coef- 
ficients necessary to extrapolate to QCD. This 
viewpoint has also been stressed in Ref. Q| in the 
context of determining the up quark mass. The 
resulting accuracy is limited by the order to which 
one calculates in ChPT. 



'Supported by DOE contract DE-FG03-96ER40956. 
t Speaker 



■^This result is also implicit in the work of Ref. j^. 
*Which we refer to as Gasser-Leutwyler (GL) coefficients. 



An important question is whether it is neces- 
sary in such a program of calculations to con- 
sider non-degenerate dynamical quarks. For the 
quantities we consider here — pion masses and de- 
cay constants — degenerate dynamical quarks suf- 
fice at one-loop in ChPT, because the expres- 
sions contain only the sum of the dynamical 
masses^. Simulations with non-degenerate dy- 
namical quarks are necessary for more accurate 
chiral extrapolations, since mass differences enter 
in two-loop ChPT results. Such simulations are 
also attractive because they bring one closer to 
QCD, and thus reduce the dependence on ChPT. 
They may also be advantageous from a computa- 
tional point of view — since a dynamical strange 
quark is relatively cheap, one may approach more 
closely to QCD by having two lighter dynamical 
quarks than three. 

For these reasons we think it likely that simula- 
tions with non-degenerate dynamical quarks will 
be done, and thus as a first step in our program 
we have calculated the one-loop expressions for 
the masses and decay constants of the pions and 
kaons in a theory with three dynamical quarks, 
only two of which are degenerate. This general- 
izes Refs. which calculated these quantities 
with any number of degenerate dynamical quarks. 

The partially quenched approximation can be 
described in terms of a Lagrangian by adding a 
number of commuting spin- 1/2 (ghost) fields §. 
For each valence quark there is a corresponding 
ghost quark, and when their masses are equal 
the valence quark loops cancel against their ghost 
counterparts. The dynamical quarks, which do 
not have corresponding ghosts, give rise to the 
quark loops that are kept in the simulation. The 
complete mass matrix is 

where V refers to valence, S to sea or dynamical, 
and a tilde indicates ghost quantities. In our the- 
ory nivi — fnvi, and msi = ms2- For baryons 
one would need to include an additional valence 
quark and its corresponding ghost. 

The approximate SUl(3) ® SUr(3) chiral sym- 
metry of QCD generalizes to the graded group 

^This appears to hold for all observables at one- loop order. 



SUl(5|2) SUr(5|2) in the PQ theory [|. As- 
suming that this larger symmetry is broken spon- 
taneously down to its vector subgroup, one finds 
a large multiplet of pseudo-Goldstone particles 
which includes both bosons (quark-antiquark and 
ghost-antighost) and fermions (quark-antighost). 
Collecting these into a matrix, <I>, and defining 
E = exp(2i$/f) and x = 2fiM, the chiral La- 
grangian describing the low-momentum interac- 
tions of these particles is |||,^ 

C = (fV4)str (asas^) - (fV4)str (xS^ + ^x) 
+La str 

+L5 str((9East) (xSt + Ex)) 
-L6Str((xSt + Ex)') 

-L8Str(xStxEt + I]xSx) +... (1) 

This form is the same as that for QCD except for 
the replacement of traces by supertraces, and for 
the appearance of the super- 77' field $0 cx str($). 
As in standard ChPT there are unknown con- 
stants at each order: f, /i, a$ and mo at leading 
order, and the GL coefficients Li at one-loop (of 
which only those that are relevant are displayed). 
Each term in Eq. ([^) can also be multiplied by an 
arbitrary even function of <i>o. These additional 
functions do not, however, enter into our final re- 
sults and are not shown. 

The QCD chiral Lagrangian is completely con- 
tained in Eq. (|l|), and can be obtained from it 
by setting valence and dynamical quark masses 
equal Moreover, since all quark mass depen- 
dence appears explicitly in Eq. (|^), and all the 
unknown coefficients are independent of masses, 
these coefficients are the same for the partially 
quenched theory and full QCD. 

Using Eq. (Q), we calculate the masses 
and decay constants of the Goldstone par- 
ticles to one- loop in ChPT [f). The re- 
sults are functions of valence and dynam- 
ical masses: M^(myi, ™si, Jriss) and 
F{mvi,mv2',msi,ms3)- Unquenched pion-like 
properties are obtained for myi — my 2 = msi, 
whereas kaon-like properties are obtained for 
TOyi — ms2 and my 2 = "^S3- The results contain 



a simple term analytic in quark masses, e.g. 

2(L4 - 2L(i){mYi + TOy2)(2msi + rngg) 

+ {L5 - 2L8){mvi + mv2?] , (2) 

and a long and unilluminating chiral logarithm 
arising from loops. Note that the dependence of 
on {mvi+rnv2f' gives the combination £5 — 
2L8, the value of which is needed to determine 
the physical m„p ||^. 

We use the full one-loop expressions to study is- 
sues that arise in the extrapolation towards QCD. 
To do this we use a simplification noted in Ref . : 
for three dynamical quarks, the super-77' has a 
mass comparable to the physical ry' and can be 
integrated out. The resulting forms are indepen- 
dent of the parameters a$ and toq and the extra 
functions of $0 implicit in Eq. (Q). They thus 
depend only on parameters that are present in 
the QCD chiral Lagrangian. We then find values 
for f, and for the GL coefficients, such that there 
are choices for the "physical" up and strange 
quark masses, denoted m„p and vrist, which re- 
produce correctly the experimental values of 
fx, m-ir and rriK- There is a range of allowed GL 
coefficients — we select representative samples by 
the additional criterion that chiral perturbation 
theory should be converging. 

Our first study concerns the importance of in- 
cluding non-analytic terms when extrapolating to 
QCD. We produce "fake" data by calculating 
and -F (using the full one-loop forms) on a pre- 
set grid of quark masses. The data is then fit to 
the one-loop analytic form predicted by ChPT 
(e.g. Eq. (|2|)), excluding the chiral logarithms. 
This fit is then used to extrapolate to the phys- 
ical values for the ratios m-^/jK and mK/fK- 
This leaves j-nj ]k as a "prediction" , which can be 
compared with the exact value. We find that, for 
grids with masses above ~ rrist/S, this procedure 
leads to 10 — 15% errors in j-^j ^k- This is true 
also if one extrapolates solely using unquenched 
data. We conclude that the extra curvature due 
to the chiral logarithms is substantial, and needs 
to be incorporated in the extrapolation. 

Most current PQ simulations use two dynam- 
ical quarks, which can be viewed as the infinite 
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Figure 2. Effect of varying the dynamical strange 
quark mass. Quark masses are my\ = m„j,, 
mv2 = rrist, msi = 17132 = ruup and rns3 = x. 
R{x) [Mk{x)/ fK{x)]/[MK{mup)/ fK{ 

strange quark mass limit of the three flavor the- 
ory. An important issue is how large an error 
one makes by excluding the dynamical strange 
quark. ChPT cannot be used to calculate this ef- 
fect, since it cannot be applied if m^j is too heavy. 
Nevertheless, by varying rrist within the range of 
validity of ChPT, one can obtain an estimate of 
the error. Fig. ^ shows how the (normalized) ra- 
tio of rriK I Jk varies as the mass of the third dy- 
namical quark is increased from rUup (i.e. three 
degenerate dynamical quarks) to rust (the phys- 
ical point). Note that the valence quark masses 
are unchanged. We see an effect of ~ 10% — large 
enough to suggest that inclusion of the dynamical 
strange quark is important. 

More details concerning these exercises, as well 
as the one-loop results for three non-degenerate 
quarks of arbitrary mass, will be given in 
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